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TECHNICAL MEMORANDUM 


HUBBLE SPACE TELESCOPE - POINTING ERROR EFFECTS 
OF NONLINEAR BALL JOINTS 

INTRODUCTION AND OBJECTIVES 


The extreme pointing accuracy of the Hubble Space Telescope (HST) requires the consideration 
of effects not normally included when performing a Pointing Control System (PCS) analysis. For this 
reason, an effort has been made to model those mechanisms which could possibly cause very small limit 
cycles, on the order of 1 milli-arc-sec. The method chosen to mount the Scientific Instrument (SI) 
Optical Benches on-ooard the HST provide one such mechanism. 

While the primary concern was for the HST, the program used for this analysis was written in 
modular form and can, therefore, be applied to other spacecraft. As shown in Figure 1, the mathematical 
model for the Optical Bench receives only the vehicle angular body rates and the spacecraft to inertial 
frame transformation. These variables are processed and returned to the main program as reactive torques 
affecting the attitude control system performance. It is assumed for purposes of this study that the 
spacecraft is operating under drag free conditions. However, an input for external forces (such as 
atmospheric resistance) has been provided to satisfy the more general situation of a vehicle experiencing 
and reacting to various external disturbances. 
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Figure 1. Optical Bench/Spacecraft simplified diagram. 

There were three primary objectives to be accomplished by this study: 

1) Define the nonlinear dynamics associated with the Si Optical Benches mounted on fre”*. ball 

joints. 

2) Develop a computer program describing these dynamics. 

3) Combine this program with an existing HST simulation to determine the effects on vehicle 
pointing accuracy. 




EQUATIONS OF MOTION 


The optical benches were assumed to be mounted as shown in Figure 2. 



Figure 2. Optical Bench mounting arrangement. 

The variables ACG, Rj, R 2 , Ry^c 316 defined * n the body-fixed frame (V j , Vj, V 3 ) to be 



ACG v j" 


R CGV1 

ACG = 

ACG V 2 

; r vcg ■ 

R CGV2 


ACG v3 


^GVS 



R 1 VI 


R 2 V 1 

R l = 

R 1V2 

; R 2 = 

R 2V2 


R 1V3 


R 2V3 

— i 


( 1 ) 


The vectors Rj and R 2 represent the distance to the geometric center of the aft and forward cavities 

(sockets), respectively. The distance from the vehicle center of mass (CM) to each cavity center and the 
optical bench CM is depicted in Figure 3 and is expressed by equation (2) 


R p = R 2 - ACG (Forward Cavity) 


( 2 ) 



R A = Ri - ACG (Aft Cavity) 


( 2 ) 

(Cone.) 


R^q = Ry^Q - ACG (Optical Bench CM) 


Let the vehicle angular rates be represented by 


Cdy = 


W V1 

«j V 2 

L u V3 


(RAD/SEC) 


(3) 



In one computational period (specified by At sec) these rates cause the center point of each cavity 
to sweep out an arc of length 


ASp — Rp GJy At 
AS a = R A GJy At . 

Expressing equation (4) in matrix notation produces 
ASp = [-Rp) [wy] At 

AS A = l-^A 1 At • 


(4) 


(5) 
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where the tilde operation for any vector 


A = 


X 

Y 

Z 


is defined to be 


0 -Z 


[A] = 


Z 

-Y 


0 -X 
v 0 


Applying the definition [equation (7)] to equation (5) gives 



0 RpV3 " R FV2 


W V1 

ASp = 

■ r PV3 0 R FV1 


W V2 


. r FV2 * R FV1 0 


_“V3 . 



0 R AV3 - R AV2 


W V1 

as a = 

- r A v3 0 R AV1 


W V2 


b R AV2 " R AV1 0 


,"V3. 


At 


At 


( 6 ) 


(7) 


( 8 ) 


Let [EV] be the vehicle to inertial frame coordinate transformation. Then the arc lengths expressed in 
inertial space are 


ASp = [EV] [-R F ] [w v l At 


AS a = [EV] t-R A ] [w v ] At 


( 9 ) 


The sum of the external and internal forces acting on the sps<«craft are represented by 
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(» 0 ) 


F = 


F V1 
F V1 
. F V3 « 


The body fixed forces given by equation (10) produce vehicle accelerations in the inertial coordinate 
frame given by 


A, = [EV]-L (ID 

My 

where My is the weight of the vehicle expressed in kilograms. Therefore, the average inertial transla- 
tional velocity of each cavity center during a single computational period is 

A I 

V AVG = V + y At < 12 > 


and the inertial distance traversed in one time period is given by 
AD = V AVG At 


(13) 


The total displacement of each cavity center can now be expressed as 


t f /At 

C F = £ (ASp + AD)j 

i=l 

t t /At 

C A - E < AS A + AD)j 

i= 1 


(14) 


where tf is the total time period being considered. 

Now consider the inertial ball displacements Bp and B A . Let the inertial frame three-axis angular 

rates of the Optical Bench be defined by (0jl. Referring to Figure 4 the arc lengths traversed by each 
ball in one time period are defined by 
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Figure 4. Translational and rotational impedance. 


AK f = [-Lpl [/§,) At 
AK a = [-L a ] [0,] At 


(15) 


Also, let the inertial translation of the Optical Bench center-of-mass in one time period be defined by the 
symbol ADy. The total displacement of each ball can now be described as 


tf/At 

Bp = 2^ (AKp + ADy )j 
i=l 

t f /At 

B A - L < ak a + ad t>, 
i=l 


(16) 


The magnitude of the force acting on each ball due to cavity wall contact is specified by 


IFpl = Kj (!Cp - Bp| - TOLp] 
i f aI = k 2 [I c a -b a |.tol a ) 


(17) 


where 


K j = spring constant for forwar’ ball. This constant represents the restoring force which acts 
when the forward ball .tacts the cavity wall. 

Kt = spring constant for aft ball. 
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TOLp “ forwird ball dead zone. This constant represents the free space (or tolerance) between 
the forward ball and the cavity wall. 

TOL A = aft ball dead zone. 


The quantity within the brackets of equation (17) represents the magnitude of the compression (exten- 
sion is not possible) of the respective springs, i.e., this is the distance the ball “sinks” into the cavity 
wall before the generated force can restore the ball to the free space. Equation (17) acts o ily after 
cavity wall contact has been made and becomes zero after contact has been lost according to equation 
(18) 


|F F I = 0 |Cp - Bp| < TOLp 

for (18) 

|f a i = 0 |C A - fl A i < TOL a 


A simplifying assumption has been made that no tangential forces can be exerted on the balls due to wall 
contact. Therefore, the forces are always directed from each ball toward the center of the cavity. The 
unit vectors expressing the direction of the forces Fp and F A are given by 


(Cp - Bp) 

Up = — 

(Cp - Bp | 



<C A - b a> 
' C A - B a‘ 


(19) 


Each ball force (expressed in the inertial frame) can now be represented as 


Fp = IFplUp 


( 20 ) 


^ - |F a |U a 


The systevr block diagram shown in Figure 5 gives a more complete definition of this logic. 
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Figure 5. System block diagram 















































CABLE CONNECTIONS 


bach Optical Bench is connected to the spacecraft through electrical cables used for power and 
scientific data transmission. These cables will have an effect on the motion of the Optical Bench since 
they have propeitie3 which produce translational and rotational spring and damping coefficients. The 
cable forces and torques arc assumed to act through the Optical Bench center of mass as depicted in 
Figure 4 Let the mass and inertia tensor be defined by M] and 


J 


1 





J V2 

0 



( 21 ) 


Establish the x, y, z coordinate frame with origin at the Optical Bench center of macs and with 
axes initially parallel to the vehicle and inertial frames. The lever arms and Lp are (referring to 
Figs. 1 and 4) 


L A = R 1 " r vcg 

l f - r 2 - Rvcg 


( 22 , 


The cable dynamics are shown in Figure 4 as an impedance (z x . z y , z z ) in each axis connecting tne 
Optical Bench center of mass and the spacecraft. 

Consider fast the translational dynamics. A mechanical representation used for equation deriva- 
tion is shown in Figure 6 and has beer simplified (for ease of drawing) such that the attach point lever 


OPTICAL BENCH 



D 


Figure 6. Translational model. 
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ami lies only along the V3 axis. Consider the vehicle to have an angular rate (wy2) about the V2 axis 
and a translation (D) along the V j axis. Under these conditions, the motion of the attach point is given 
by the sum of D and D^y where 


D RT = / V RT dt ~ ” /^CG W V2 dt (23) 

Expressing this in the more general three-dimensional case and in matrix notation gives 

Drt = V RT = I _R CG^ w V • (24) 


Writing the transfer function from the input force F to the attach point produces, 

F = M j by + dy (Dy - Dj^y -D) + Ky (Dy - Dpy - D) . (25) 

Rearranging equation (25) and introducing LaPlace notation gives 

F = (MjS“ + dyS + Ky)Dy - (dyS + Ky) Dj^y - (dyS + Ky) D . (26) 

If it is assumed that the vaiiibles D^y and D are zero (as would be the case for a stationary vehicle) 
then equation (26) simplifies to 

Dy 1 /M j 1 /M i 

-r = - ? --5 t • ( 27 > 

t S 2 + dy/Mj S + Ky/Mj S 2 + 2f o; N S + W N 2 

For the purpose of this report it was assumed that the spring constant (Ky) and damping ratio (J) were 

Ky = 1 Ib/in. = 175.1 18 N/m and J = 0.005 (unitless). 

With these assumptions the natural frequency and damping coefficient are calculated to be 
f = 0. 1 74 Hz and dy=1.6N-S/m 

Since the Optical Bench mass is invariant all of the translational coefficients will be the same for all 
vehicle axes 

A simplified vie 'f the n*. hanical model used for rotational equation derivation is shown in 
Figure 7. 
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V 3 



Figure 7. Rotational model. 

The torque applied to the Optical Bench inertia [equation (21)] due to wall contact is 

T= [Lp] [F F ] + [L a ] [F a ] , (28) 

where and Lp are defined in Figure 4. 

Summation of the torques on each inertia produces the following equation, 


T = Jj &Z + dR (Pz ~ + (fi - 9) 

where 



V 


Q v\ 

01 = 

0V2 

; e = 

d V2 


0V3 


d V2 


are the rotation angles of the Optical Bench and the vehicle, respectively. 

Again introducing LaPlace notation and letting u>\ = 6 = 0 produces 


(29) 


(30) 
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( 31 ) 


0 

T 



S 2 + d R /Jj S + Kr/Jj 



+ 2ftu^S + 


With the assumption that 


K R 


1.3558 


Nm 

RAD 


f = 0.005 unitless 


the following constants are calculated 



“0.0877" 



“0.0246 

II 

cc 

T3 

0.1431 

(NmS) ; 

f N “ 

0.015! 


0.1431 



0.0151 


(Hz) 


The forces and Fp are calculated in the inertial frame; however, the translational and rota- 
tional spring and damping constants only have meaning in the vehicle reference axis system. Conse- 
quently, the forces and torques generated as a result of motional differences between the vehicle and 
Optical Bench must be converted into the vehicle frame and back to the inertial frame. Equation (25) 
can be written as 


^ IF + F DAMP + F S1 


(32) 


where 


F DAMP ~ d T ( ^ + ^RT ‘ Dj) " d T< V “ V T + V RT> 


(33) 


Fg — Ky(D — Dj + Dr"p) 


Since the spring and damping factors can be applied only in the vehicle frame, equation (33) must be 
rewritten as 


12 


f DAM? * l EV 1 {td T l x UvE] IV - V T 1 + V RT )} 

F S 1 IEV1 | (K X J x ([VE] ID - Dj.) + D RT )J 

where (VE) = [EV]^ = vehick to inertial frame transformation. 
For the rotational dynamics, equation (29) is expressed as 

0 = J l _1 lT + T DAMP + T sl 

where 

T DAMP = d Rl w V * fa 
T s = K R [0 - (31 


(36) 


Again converting to the vehicle frame and back to the inertial reference gives 


t DAMP ■ [EV] (d R x <«v - IVE] [0,1 )} 

T S =[EV| |k r x(6 v -[VE] [?,]) l . 

Equations (34) and (37) represent the dynamics as implemented and as shown in the block diagram 
(Fig. 5). 

Referring to equation (11) the vehicle acceleration can now be expressed in more detail by 


(37) 


(EV) [F] [^a + ^F + ^DAMP + Eg " (EXTERNAL FORCES)] 

^ My My 

It is possibk for the Optical Bench to transmit reactive torques to the vehicle through three 
different paths: 

1) Forward and aft wall contact forces acting through their respective lever arms. 

2) Optical Bench rotational damping and spring torques. 

3) Translational spring and damping forces acting at the Optical Bench attach point. 
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The total reactive torque (in vehicle coordinates) can now be expressed as 


T r = -[VE] | [Rp] [F p ] + [R a ] 1F a ] + [Reel l p DAMP + F S^ * T DAMP + T S } • < 39 > 

Equation (39) represents the reactive torques that drive the spacecraft attitude control system. 


SIMPLIFYING ASSUMPTIONS 


All the necessary equations have been derived and the loop can now be closed through equation 
(39) as shown in Figure 5. Certain simplifying assumptions have been made in the equation derivation 
either because their effects were considered to be negligible or because the data for their implementation 
was not available. 

1) The arc lengths AS, ADy, and AK were considered to be s^aight lines instead of curved paths. 

2) No tangential forces are exerted on the Optical Bench balls by cavity wall contact. Contact 
forces are always directed along a line through the contact point and the cavity center. 

3) Cable damping and spring constants for rotation and translation act through the Optical Bench 
center of mass. 

4) The Optical Bench inertia tensor was considered to be diagonal (i.e., all coupling terms were 

zero). 


5) The ball and socket dead zones (specified by TOLp and TOL A ) are small enough to allow 
the Optical Bench reference frame to ■'e considered congruent with the vehicle frame. 

6) The total torque on the Optical Bench is expressed by 


T t = Jj fi + (fi x H) 


where H is the angular momentum. Assuming a diagonal inertia tensor and small angular rates, the term 
(0 x H) can be ignored. 

7) The ball contact spring stiffness was assumed to be constant. 

8) The Optical Bench was assumed to be rigid. 
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INPUT DATA AND PROGRAM RESULTS 


The Faint Object Camera (FOC) is the only Scientific Instrument on beard the HST which has 
the ball and socket mounting arrangement at each end. Referring to Figure 2 the variables shown for 
the FOC are given by (all units are in meters) 



6.617 


4.461 


3.1 


5.125 

ACG = 

0.005 

: R VCG " 

-0.297 

; R,= 

-0.102 

; R 2 = 

-0.7 


-0.145 


-0.325 


-0.102 


-0.7 


Combining the data in Figure 2 according to previously developed equations produces the variables of 
Figure 3 



-3.517 


-2.156 


-1.492 

r a = 

-0.i07 

I ^G ' 

-0.302 

; Rp = 

-0.705 


0.043 


-0.18 


-0.555 


and Figure 4 



~1.36T 


"~0.66? 

n 

< 

-j 

0.195 

ii 

u. 

-0.403 


0.233 


-0.375 


The FOC mass and inertia terms are defined by 


and 


M, = 146 Kg 


J l = 


7.256 0 0 

0 39.207 0 


0 


0 


Kg m' 


39.207 


where J j is given about the FOC center of mass. 


As a test case the vrhicle was maneuvered 0.05 deg about the positive V 2 axis. This rotation was 

sufficient to cause an aft cavity wall collision. The system performance was examined with and without 
flexible body dynamics. The rotation about the V 2 (maneuver axis) produced very similar results in each 

case as shown in Figures 8 and 9. Each plot shows a rotation to 0.05 deg (180 aro-sec) but the scale is 
such that the small variations due to body bending are not evident. Plotting the angular error in the 
maneuver axis produces Figures 10 and 1 1. The maneuver takes approximately 80 sec to complete. 

The small variations in each plot are less than 1 milli-arc-sec and are principally due to the translational 
and rotational dynamics.of the suspended mass. Flexible body effects are much more pronounced in 
Figures 1 2 through 1 5 which represent the angular rotations for the V j and axes. Collisions between 

the aft ball and the cavity wall are most evident in Figures 14 and 15 which indicate six collisions of the 
aft ball at 14, 22, 26, 54, 62, and 67 sec. After the maneuver has been completed, the flexible body 
effects are once again very small and the variations can be attributed primarily to mass translation 


CONCLUSIONS 


The simulation results indicate that there is no strong coupling from the wall/ball collisions into 
the vehicle elastic body. Flexible body “ringing” soon damps out after the vehicle rotation has been accomp- 
lished and the last collision has occurred. Pointing axis deviations as determined from the RMS value of the V 2 

and V 3 axis rotations are less than 1 x 10 arc-sec. The principal cause of the damped vibration seen 

beginning at 80 sec is the translational effects of the Optical Bench acting through the cable connections. 
Studies will continue to determine if uncertainties in flexible body and translational mode frequencies 
could cause an instability or limit cycle condition to exist. 
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Figure 10. V 2 axis error with rigid body dynamics (difference between 
commanded and actual vehicle angle). 
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Figure 11. V 2 axis error with flexible body dynamics (difference between 
commanded and actual vehicle angle). 
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